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Abstract 

We study the production and decays of top squarks (stops) at the Tevatron collider in models of 
low-energy supersymmetry breaking. We consider the case where the lightest Standard Model (SM) super- 
partner is a light neutralino that predominantly decays into a photon and a light gravitino. Considering 
the lighter stop to be the next-to-lightest Standard Model superpartner, we analyze stop signatures associ- 
ated with jets, photons and missing energy, which lead to signals naturally larger than the associated SM 
backgrounds. We consider both 2-body and 3-body decays of the top squarks and show that the reach of 
the Tevatron can be significantly larger than that expected within either the standard supergravity models 
or models of low-energy supersymmetry breaking in which the stop is the lightest SM superpartner. For a 
modest projection of the final Tevatron luminosity, C ~ 4 fb _1 , stop masses of order 300 GeV are accessible 
at the Tevatron collider in both 2-body and 3-body decay modes. We also consider the production and 
decay of ten degenerate squarks that are the supersymmetric partners of the five light quarks. In this case 
we find that common squark masses up to 360 GeV are easily accessible at the Tevatron collider, and that 
the reach increases further if the gluino is light. 
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1 Introduction 



The Standard Model with a light Higgs boson provides a very good description of all experimental data. 
The consistency of the precision electroweak data with the predictions of the Standard Model suggests 
that, if new physics is present at the weak scale, it is most probably weakly interacting and consistent with 
the presence of a light Higgs boson in the spectrum. Extensions of the Standard Model based on softly 
broken low energy supersymmetry (SUSY) provide the most attractive scenarios of physics beyond 
the Standard Model fulfilling these properties. If the supersymmetry breaking masses are of the order of 
1 TeV, supersymmetry stabilizes the hierarchy between the Planck scale Mp and the electroweak scale. 
Furthermore, the minimal supersymmetric extension of the Standard Model (MSSM) significantly improves 
the precision with which the three gauge couplings unify and leads to the presence of a light Higgs boson 
with a mass below 135 GeV 

Perhaps the most intriguing property of supersymmetry is that local supersymmetry naturally leads 
to the presence of gravity (supergravity) . In the case of global supersymmetry, a massless spin one-half 
particle, the Goldstino, appears in the spectrum when supersymmetry is spontaneously broken. In the case 
of local supersymmetry, the Goldstino provides the additional degrees of freedom to make the gravitino a 
massive particle ||. In the simplest scenarios, the gravitino mass tHq is directly proportional to the square 
of the supersymmetry breaking scale V -^susy : 

ro<3 ~ F SVSY /V3M P . (1.1) 

The relation between the supersymmetry breaking scale V -Fsusy and the supersymmetry particle 
masses depends on the specific supersymmetry breaking mechanism. In general, the superpartner masses 
-^SUSY are directly proportional to -Fsusy and inversely proportional to the messenger scale M m at which 
the supersymmetry breaking is communicated to the visible sector: 

Msusy * (1.2) 

where a is the characteristic coupling constant coupling the messenger sector to the visible one. If the 
breakdown of supersymmetry is related to gravity effects, M m is naturally of the order of the Planck scale 
and q is of order one, so that Mstjsy is at the TeV scale if V -^susy ~ 10 11 GeV. In gauge mediated 
scenarios (GMSB) H,^], instead, the couplings a are associated with the Standard Model gauge couplings 
times a loop suppression factor, so that values of Fs\jsy/M m of order 100 TeV yield masses for the lighter 
Standard Model superpartners of order 100 GeV. 

When relevant at low energies, the gravitino interactions with matter are well described through the 
interactions of its spin 1/2 Goldstino component ||]. The Goldstino has derivative couplings with the 
visible sector with a strength proportional to l/fsuSY- In scenarios with a high messenger scale, of order 
Mp, Eqs. |0| and [E^ imply that the gravitino has a mass of the same order as the other SUSY particles, 
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and its interactions are extremely weak. In such scenarios, the gravitino plays no role in the low-energy 
phenomenology. However, in low energy supersymmetry breaking scenarios such as GMSB in which the 
messenger scale is significantly lower than the Planck scale, the supersymmetry breaking scale is much 
smaller. Typical values in the GMSB case are M m ~ 10 5 — 10 8 GeV, leading to a supersymmetry breaking 
scale V -P31TSY roughly between 10 5 and a few times 10 6 GeV. The gravitino then becomes significantly 
lighter than the superpartners of the quarks, leptons and gauge bosons, and its interaction strength is 
larger. As the lightest supersymmetric particle (LSP), the gravitino must ultimately be produced at the 
end of all superparticle decay chains if -R-parity is conserved (for an analysis of the case of i?-parity violation 
see Ref. [§). 

Depending on the strength of the gravitino coupling, the decay length of the next-to-lightest supersym- 
metric particle (NLSP) can be large (so that the NLSP is effectively stable from the point of view of collider 



phenomenology; this occurs when V -Fsusy <^ 1000 TeV), intermediate (so that the NLSP decays within the 
detector giving rise to spectacular displaced vertex signals; this occurs when 1000 TeV > yj Fsusy ^ 100 



TeV), or microscopic (so that the NLSP decays promptly; this occurs when V Fstjsy ^ 100 TeV) The 
decay branching fractions of the Standard Model superpartners other than the NLSP into the gravitino 



are typically negligible. However, if the supersymmetry breaking scale is very low, \J Fsusy <C 100 TeV 
(corresponding to a gravitino mass <C 1 eV), then the gravitino coupling strength can become large enough 
that superpartners other than the NLSP may decay directly into final states containing a gravitino H 10 1 . 
In any case, the SUSY-breaking scale must be larger than the mass of the heaviest superparticle; an 
approximate lower bound of \J -Fsusy — 1 TeV yields a gravitino mass of about 10" 3 eV. 

In many models, the lightest supersymmetric partner of a Standard Model particle is a neutralino, Xi- 
The decay widths of the Xi m ^>o various final state particles are given by: 

4 / \ ^ 

r« S XG) , KxNx ^ - ^) . 

where Kx is a projection factor equal to the square of the component in the NLSP of the superpartner of 
the particle X, and Nx is the number of degrees of freedom of X. If X is a photon, for instance, Nx = 2 
and 

K x = \N u cos9 w + N 12 sm9 w \ 2 , (1.4) 

where Nij is the mixing matrix connecting the neutralino mass eigenstates to the weak eigenstates in the 
basis B,W,H ly H 2 . 

If the neutralino has a significant photino component, it will lead to observable decays into photon and 
gravitino. Since the heavier supersymmetric particles decay into the NLSP, which subsequently decays 
into photon and gravitino, supersymmetric particle production will be characterized by events containing 
photons and missing energy. This is in contrast to supergravity scenarios, where, unless very specific 
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conditions are fulfilled [11, 12| , photons do not represent a characteristic signature. The presence of two 
energetic photons plus missing transverse energy provides a distinctive SUSY signature with very little 
Standard Model background. 

If Xi decays with a large branching ratio into photons, its mass is severely constrained by LEP data. 
However, this bound is model-dependent: there is no tree- level coupling of the 7 and Z to two binos or 
to two neutral winos, but LEP could produce bino or neutral wino pairs through i-channel exchange of 
selectrons. Since the bino is associated with the smallest of all gauge interactions and, in addition, its mass 
is more strongly renormalized downward at smaller scales compared to the wino mass, in many models the 
lightest neutralino has a significant bino component. Therefore, if the NLSP is approximately a pure bino, 
the bounds on its mass depend strongly on the selectron mass. For selectron masses below 200 GeV, the 
present bound on such a neutralino is of about 90 GeV |13| . Then, as emphasized before, it will decay via 
Xi - > jG. If Xi is heavy enough, then the decays Xi ~^ %G and Xi ~^ h°G are also kinematically allowed; 
however, the decay widths into these final states are kinematically suppressed compared to the 7G final 
state and will only be important if the photino component of Xi is small or if Xi is significantly heavier 
than Z and h° @. 

In most SUSY models it is natural for the lighter top squark, t\, to be light compared to the other 
squarks. In general, due to the large top Yukawa coupling there is a large mixing between the weak 
eigenstates II and £r, which leads to a large splitting between the two stop mass eigenstates. In addition, 
even if all squarks have a common mass at the messenger scale, the large top Yukawa coupling typically 
results in the stop masses being driven (under renormalization group evolution) to smaller values at the 
weak scale. An extra motivation to consider light third generation squarks comes from the fact that light 
stops, with masses of about or smaller than the top quark mass, are demanded for the realization of the 
mechanism of electroweak baryogenesis within the context of the MSSM [ |i~4"| . 

In this paper, we examine in detail the production and decay of top squarks at Run II of the Tevatron 
collider in low-energy SUSY breaking scenarios, by assuming that the NLSP is the lightest neutralino, 
which decays promptly to jG. We also investigate the production and decay of the other squarks, and 
provide an estimate of the reach of Run II of the Tevatron in the heavy gluino limit. We work in the context 
of a general SUSY model in which the SUSY particle masses are not constrained by the relations predicted 
in the minimal GMSB models. We assume throughout that the gravitino coupling is strong enough (or, 
equivalently, that the scale of SUSY breaking is low enough) that the NLSP decays promptly. This implies 
an upper bound on the supersymmetry breaking scale of a few tens to a few hundred TeV depending 
on the mass of the NLSP. Our analysis can be extended to higher supersymmetry breaking scales for which 
the NLSP has a finite decay length, although in this case some signal will be lost due to NLSP decays 
outside the detector. At least 50% of the diphoton signal cross section remains though for NLSP decay 
lengths ct < 40 cm M; this corresponds to a supersymmetry breaking scale below a few hundred to about 
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a thousand TeV, depending on the mass of the NLSP. Moreover, the displaced vertex associated with a 
finite decay length could be a very good additional discriminator for the signal. Thus, in totality, our 
choice is certainly not an overly optimistic one. 

This paper is organized as follows. In Sec. |2| we review the stop pair production cross section at the 
Tevatron. In Sec. || we summarize previous studies of stop production and decay at Run II of the Tevatron. 
In Sec. |] we discuss the SUSY parameter space and the relative partial widths of the various stop decay 
modes. In Sec. |5] we describe the signal for each of the stop decay modes considered. We describe the 
backgrounds and the cuts used to separate signal from background in each case, and give signal cross 
sections after cuts. This gives the reach at Run II. We also note the possibility that stops can be produced 
in the decays of top quarks. In Sec. ^ we consider the production and decay of 10 degenerate squarks that 
are the supersymmetric partners of the five light quarks. In Sec. ^ we summarize our conclusions. 



2 Top squark production at Tevatron Run II 

Top squarks are produced at hadron colliders overwhelmingly via the strong interaction, so that the tree- 
level cross sections are model independent and depend only on the stop mass. The production modes of 
the lighter stop, t±, at the Tevatron are qq — ► t\t\ and gg — > tit*. The cross sections for these processes 
are well known at leading order (LO) p£, 16|, and the next-to-leading order (NLO) QCD and SUSY-QCD 



corrections have been computed |17J and significantly reduce the renormalization scale dependence. The 
NLO cross section is implemented numerically in PROSPINO [17, 18|. 



We generate stop events using the LO cross section evaluated at the scale n = rrif, improved by the 
NLO K-factor [] obtained from PROSPINO |17,lq| (see Fig. |l|). The K-factor varies between 1 and 1.5 for 



mi decreasing from 450 to 100 GeV. We use the CTEQ5 parton distribution functions [19|. We assume 
that the gluino and the other squarks are heavy enough that they do not affect the NLO cross section. 
This is already the case for gluino and squark masses above about 200 GeV |17[| . 

Top squarks can also be produced via cascade decays of heavier supersymmetric particles, with a highly 
model dependent rate. To be conservative, we assume that the masses of the heavier supersymmetric 
particles are large enough that their production rate at Tevatron energies can be neglected. 



3 Previous studies of stops at Run II 

A number of previous studies have considered the prospects for stop discovery at Run II of the Tevatron, 
which we summarize here. In general, the most detailed SUSY studies have been done in the context of 
1 Although gluon radiation at NLO leads to a small shift in the stop pr distribution to lower pr values fll7| , |l8| , we do not 
expect this shift to affect our analysis in any significant way. 
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Figure 1: LO and NLO cross sections for stop pair production in pp collisions at Tevatron Run I (1.8 TeV) 



and Run II (2.0 TeV), from PROSPINO |jT7 , 18|. Cross sections are evaluated at the scale /i = m^. 
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supergravity; in this case SUSY is broken at the Planck scale so that the gravitino plays no role in the 
collider phenomenology. Then the lightest neutralino is the LSP and ends all superparticle decay chains. 
The reach of the Tevatron for a number of stop decay modes has been analyzed in Refs. [ 2Cj , 21 1. The signal 
depends on the decay chain, which in turn depends on the relative masses of various SUSY particles. For 
sufficiently heavy stops, the decay t — » tx\ will dominate. This channel is of limited use at Run II because 
the stop pair production cross section falls rapidly with increasing stop mass, and this channel requires 
m ( - > mj + m x l l ' which is quite heavy for the Tevatron in the case of minimal supergravity.^ For lighter 
stops, if a chargino is lighter than the stop then t — ► bxf tends to dominate (followed by the decay of the 
chargino). The details of the signal depend on the Higgsino content of xt • ^ Xi has a mass larger than 
mi — nib, the previous decay does not occur and the three-body decay i — > bW + Xi dominates; this decay 
proceeds through the exchange of a virtual top quark, chargino, or bottom squark. If the stop is too light 
to decay into an on-shell W boson and Xi> then the flavor-changing decay t — > ex? tends to dominate. 
Finally, if a sneutrino or slepton is light, then t — > bl + vi or t —* b£ + V£, respectively, will occur (followed 
by the decays of the slepton or the sneutrino, if it is not the LSP). At Run II with 2 (20) fb _1 of total 
integrated luminosity, in the context of minimal supergravity one can probe stop masses up to 160 (200) 
GeV in the case of the flavor changing decay, while stop masses as high as 185 (260) GeV can be probed 



if the stop decays into a bottom quark and a chargino [21|. A similar reach holds in the case of a light 
sneutrino |^l| and in the case of large tan/3 when the stop can decay into bruxi [22|. 



One can also search for stops in the decay products of other SUSY particles |2l[]. In top decays, the 
process t — ► tg is already excluded because of the existing lower bound on the gluino mass 0. Other 
possibilities are x~ ~* bt* and g —* tt*. Finally, the decays b — > tW~ and tH~ have to compete with the 
preferred decay, b — > bxl, and so may have small branching ratios (depending on the masses of t, b, X\ an d 
H~). The signals for these processes at the Tevatron Run II have been considered in minimal supergravity 
in Ref. ||. 

Relatively few studies have been done in the context of low-energy SUSY breaking with a gravitino 
LSP. A study of GMSB signals performed as part of the Tevatron Run II workshop @ considered the 
decays of various SUSY particles as the NLSP. As discussed before, the NLSP in such models will decay 
directly to the gravitino and Standard Model particles. If the stop is the NLSP in such a model, then it 
will decay via t — * t^*'G — > bW + G (for > + mw)- Note that because G is typically very light in 



such models, > mw + mb is sufficient for this decay to proceed with an on-shell W boson. Ref. [21| 
found sensitivity at Run II to this decay mode for stop masses up to 180 GeV with 4 fb _1 . This stop 
decay looks very much like a top quark decay; nevertheless, even for stop masses near mt, such stop decays 
2 In low-energy SUSY breaking scenarios, however, the mass range > m t + m^o is interesting at Tevatron energies 

because the distinctive signal allows backgrounds to be reduced to a very low level, as we will show. 

3 The bound on the gluino mass is, however, model dependent. Under certain conditions, an allowed window exists for 

gluino masses below the gauge boson masses |23lp 
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can be separated from top quark decays at the Tevatron using kinematic correlations among the decay 
products |25|j . 

Finally, Ref. || considered general GMSB signals at the Tevatron of the form 77 $t + X. The authors 
of Ref. H provide an analysis of the possible bounds on the stop mass coming from the Run I Tevatron 
data. They analyze the stop decay mode into a charm quark and a neutralino, and also possible three 
body decays, by scanning over a sample of models. They conclude that stop masses smaller than 140 GeV 
can be excluded already by the Run I Tevatron data within low energy supersymmetry breaking models 
independent of the stop decay mode, assuming that m^o > 70 GeV. 

4 Top squark decay branching ratios 

The decay properties of the lighter stop depend on the supersymmetric particle spectrum. Of particular 
relevance are the mass splittings between the stop and the lightest chargino, neutralino and bottom squark. 
In our analysis, we assume that the charginos and bottom squarks are heavier than the lighter stop, so 
that the on-shell decays t — > bW and t — > x + 6 are kinematically forbidden. Then the details of the stop 
decay depend on the mass splitting between the stop and the lightest neutralino. If < m\y + mt, + m^o, 
two decay modes are kinematically accessible: 

1. the flavor-changing (FC) two-body decay t — > cx\- This two-body decay proceeds through a flavor- 
changing loop involving W + , H + or x + exchange or through a tree- level diagram with at — c mixing 
mass insertion; 

2. the four-body decay via a virtual W boson, t — ► W + *bxi — ► jjbxi or ivb\\ p6fl - 

The branching ratio of the stop decay into charm and neutralino strongly depends on the details of the 
supersymmetry breaking mechanism. In models with no flavor violation at the messenger scale, the whole 
effect is induced by loop effects and receives a logarithmic enhancement which becomes more relevant for 
larger values of the messenger mass. In the minimal supergravity case, the two-body FC decay branching 
ratio tends to be the dominant one. In the case of low energy supersymmetry breaking, this is not necessarily 
the case. Since the analysis of the four-body decay process is very similar to the three-body decay described 
below for larger mass splittings between the stop and the lighter neutralino, here we shall analyze only the 
case in which the two-body FC decay t — ► cx± is the dominant one whenever mi < mw + fnb + m x1 • 

For larger mass splittings, so that mw + rrib + m^o < mi < m t + m^o, the three-body decay t — ► W + b\\ 
becomes accessible, with \i ~^ lG. This stop decay proceeds through a virtual top quark, virtual charginos, 
or virtual sbottoms. Quite generally, this 3-body decay will dominate over the 2-body FC decay in this 
region of phase space. 
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For still heavier stops, rtif > mt + m^o, ^ ne two-body tree-level decay mode t — ► tx\ becomes kinemat- 
ically accessible and will dominate. Although the 3-body and 2-body FC decays are still present, their 
branching ratios are strongly suppressed. 

Let us emphasize that, since the bino is an admixture of the zino and the photino, a pure bino neutralino 
can decay into either 7G or ZG (see Eq. |1.3| ). If the lightest neutralino is a mixture of bino and wino 
components, then the relative zino and photino components can be varied arbitrarily, leading to a change 
in the relative branching ratios to 7G and ZG. If the lightest neutralino contains a Higgsino component, 
then the decay to h°G is also allowed. We show in Fig. [2] the branching ratio of the lightest neutralino into 
7G as a function of its mass and Higgsino content. 



5 Top squark signals in low-energy SUSY breaking 

As discussed in the previous section, the decay properties of the lighter stop depend primarily on the mass 
splitting between the stop and the lightest neutralino. In this section, we proceed with the phenomeno- 
logical analysis of the signatures of top squark production associated with the different decay modes. In 



Sec. 5.1 we shall analyze the signatures associated with the two-body FC decay, which after the neutralino 



decay leads to t — > C7G. In Sec. 5.2 we shall analyze the signatures associated with the three-body de- 
cay, which after neutralino decay leads to t — ► bW + jG. The two-body decay t — > txi, which typically 
dominates for > nit + m x°i leads to the same final state as the three-body decay and therefore the 



discussion of this case will be included in Sec. [5^. Finally, in Sec. 5.3 we consider the possibility that stops 
are produced in the decays of top quarks. 

5.1 Two-body FC decay: t cyG 

With the stop undergoing the aforementioned decay, the final state consists of a pair each of charm jets, 
photons and (invisible) gravitinos. As we will show below, the backgrounds to this process are small enough 
that we need not require charm tagging. Thus the signal consists of 

2(jets) + 2 7 + pr. 

The selection criteria we adopt are: 

1. each event must contain two jets and two photons, each of which should have a minimum transverse 
momentum (j?t > 20 GeV) and be contained in the pseudorapidity range —2.5 < n < 2.5; 

2. the jets and the photons should be well separated from each other; namely, 

SR j:j > 0.7, <5i? 77 > 0.3, 8R h > 0.5 
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Figure 2: Branching ratio for the decay of the lightest neutralino into a photon and a gravitino as a 
function of the neutralino mass. Shown are the branching ratio if the neutralino is a pure bino (solid 
line) and for 20% and 50% Higgsino admixtures (long and short dashed lines, respectively) assuming that 
rrtyfi = 120 GeV and that the other MSSM Higgs bosons are very heavy. For the Higgsino admixture in 
the lightest neutralino, we choose the H1-H2 mixing so that the field content is aligned with that of h° 
and the longitudinal component of the Z boson in order to minimize the branching ratio to photons. 
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Background 


Cross section after cuts 


after 7 ID 


jjllZ, Z -» vv 


~ 0.2 fb 


~ 0.13 fb 


jj~fvv + 7 radiation 


~ 0.002 fb 


~ 0.001 fb 


6677, CC77 


< 0.1 fb 


~ 0.06 fb 


nil 


~ 0.2 fb 


~ 0.13 fb 


Backgrounds with fake photons: 






jj(ee -» 77) 


~ 5 x 10 fb 


~ 5 x 10 fb 


33l(j -> 7) 


~ 0.8 fb 


~ 0.8 fb 


33(33 -> ll) 


~ 0.8 fb 


~ 0.8 fb 


Total 


~ 2 fb 


~ 2 fb 



Table 1: Backgrounds to tt* — > jj'77 ^t- The photon identification efficiency is taken to be e 7 
each real photon. See text for details. 



0.8 for 



where 5R 2 = 5r] 2 + <5</> 2 , with 5rj (5(f)) denoting the difference in pseudorapidity (azimuthal angle) of 
the two entities under consideration; 

3. the invariant mass of the two jets should be sufficiently far away from the W- and Z-masses: 
m-jj £ (75 GeV, 95 GeV); 

4. each event should be associated with a minimum missing transverse momentum (j/t > 30 GeV). 

The photons and jets in signal events tend to be very central; in particular, reducing the pseudorapidity 
cut for the two photons to —2.0 < r/ < 2.0 would reduce the signal by less than about 3%. (This change 
would reduce the background by a somewhat larger fraction.) Apart from ensuring observability, these 
selection criteria also serve to eliminate most of the backgrounds, which are listed in Table |l[ 

A primary source of background is the SM production of jjn^i^i where the jets could have arisen from 
either quarks or gluons in the final state of partonic subprocesses. A full diagrammatic calculation would 
be very computer-intensive and is beyond the scope of this work. Instead, we consider the subprocesses 
that are expected to contribute the bulk of this particular background, namely pp — ► 2j + 27 + Z + X with 
the Z subsequently decaying into neutrinos. These processes are quite tractable and were calculated with 



the aid of the helicity amplitude program Madgraph \27\. On imposition of the abovementioned set of 
cuts, this background at the Run II falls to below ~ 0.2 fb. 

An independent estimate of the jjn^iVi background may be obtained through the consideration of the 
single-photon variant, namely jjlViVi production, a process that MADGRAPH can handle. After imposing 
the same kinematic cuts (other than requiring only one photon) as above, this process leads to a cross 
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section of roughly 0.2 fb. Since the emission of a second hard photon should cost us a further power of 
a em , the electromagnetic coupling constant, this background falls to innocuous levels. We include both 
this estimate and the one based on Z production from the previous paragraph in Table |l|. While a naive 
addition of both runs the danger of overcounting, this is hardly of any importance given the overwhelming 
dominance of one. 

A second source of background is 6677 (0x77), with missing transverse energy coming from the semilep- 
tonic decay of one or both of the b (c) mesons. In this case, though, the neutrinos tend to be soft due 
to the smallness of the b and c masses. Consequently, the cut on $t serves to eliminate most of this 
background, leaving behind less than 0.1 fb. This could be further reduced (to <J 0.001 fb) by vetoing 
events with leptons in association with jets. However, such a lepton veto would significantly impact the 
selection efficiency of our signal, which contains cc, thereby reducing our overall sensitivity in this channel. 

A third source of background is jj'77 production, in which the jet (light quark or gluon) and/or photon 
energies are mismeasured leading to a fake ftp. To simulate the effect of experimental resolution, we 



use a (very pessimistic) Gaussian smearing: 5Ej/Ej = 0.1 + 0.6/ J Ej(GeV) for the jets and 5E^/E^ 



0.05 + 0.3/y Ey(GeV) for the photons.^ While the production cross section is much higher than that of 
any of the other backgrounds considered, the ensuing missing momentum tends to be small; in particular, 
our cut on Et reduces this background by almost 99% |] On imposition of our cuts, this fake background 
is reduced to ~ 0.2 fb.f] 

Finally, we consider the instrumental backgrounds from electrons or jets misidentified as photons. 



Based on Run I analyses [29] of electron pair production and taking the probability for an electron to fake 
a photon to be about 0.4%, we estimate the background due to electrons faking photons to be of order 
5 x 10~ 4 fb. More important is the background in which a jet fakes a photon. Based on a Run I analysis |28| 
and taking the probability for a jet to fake a photon to be about 0.1%, we estimate the background due 
to jjj'7 in which one of the jets fakes a second photon to be about 0.8 fb. Similarly, we estimate that the 
background due to jjjj in which two of the jets fake photons is somewhat smaller; to be conservative we 
take it to be of the same order, i.e., 0.8 fb. 

Having established that the backgrounds are small, let us now turn to the signal cross section that 
survives the cuts. In Fig. |2| we present these as contours in the m^-m^o plane. We assume that the 
branching ratio of t — ► ex® dominates in the region of parameter space that we consider here. The 
4 We have also performed similar smearing for the other background channels as well as for the signal. However, there it 

hardly is of any importance as far as the estimation of the total cross section is concerned. 

5 This reduction factor is in rough agreement with that found for 7 + j events in CDF Run I data in Ref. [gSJ . 

6 We note that in the squark searches in supergravity scenarios the signal is jj $t\ in this case a similar background due 

to dijet production with fake j/t is present. This background is larger by two powers of a than the .7,777 background, yet it 

can still be reduced to an acceptable level by a relatively hard cut on j/t (see, e.g., Ref. El). 
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branching ratio of Xi ~^ lG is taken from Fig. [2| assuming that Xi is a pure bino.[] All our plots are 
made before detector efficiencies are applied. With a real detector, each photon is identified with about 
e 7 = 80% efficiency; thus the cross sections shown in Fig. |3] must be multiplied by = 0.64 in order to 
obtain numbers of events.^] While the production cross sections are independent of the neutralino mass, 
the decay kinematics have a strong dependence on m^o. For a given mj, a small mass splitting between 
?7ij and m^o would imply a soft charm jet, which would often fail to satisfy our selection criteria. This 
causes the gap between the cross section contours and the upper edge of the parameter space band that we 
are exploring in Fig. ^. This is further compounded at large m^o by the fact that a large neutralino mass 
typically implies a smaller Xi ~~ lG branching fraction (see Fig|2|). On the other hand, a small m^o results 
in reduced momenta for the gravitino and the photon, once again resulting in a loss of signal; however, 
this is important only for m^o <> 70 GeV, which is not relevant in this search channel. 

The signal cross sections are fairly substantial. In particular, the dark area in Fig. ^ corresponds to 
a signal cross section of 50 fb or larger at Run I of the Tevatron. Taking into account the identification 
efficiency of 64% for the two photons, such a cross section would have yielded 3 signal events in the 100 
pb _1 collected in Run I over a background of much less than 1 event. Run I can thus exclude this region 
at 95% confidence level. In particular, we conclude that Run I data excludes stop masses up to about 
200 GeV, for large enough mass splitting between the stop and the neutralino. When the stop-neutralino 
mass splitting is small (i.e., less than about 10 — 20 GeV), the charm quark jets become too soft and the 



signal efficiency decreases dramatically. For comparison, a DO search [ 30 1 for inclusive pp —* x.2 + X with 

X2 - ► 7Xi m the context of minimal supergravity yields a limit on the production cross section of about 1 

pb for parent squark masses of order 150-200 GeV. Interpreting this in terms of stop pair production with 

X2 — * 7Xi reidentified as Xi ~~ * lG increases the signal efficiency by a factor of ~ 2.7 because every SUSY 

event now contains two photons the DO analysis then yields a stop mass bound of about 180 GeV, in 

rough agreement with our result .0 In addition, Ref. 0j projected a Run I exclusion of stops in this decay 

channel for masses below about 160-170 GeV, again in rough agreement with our result. 

To claim a discovery at the 5a level, one must observe a large enough number of events that the 

probability for the background to fluctuate up to that level is less than 5.7 x 10 -7 . Because the number of 

expected background events in this analysis is small, we use Poisson statistics to find the number of signal 

events required for a 5er discovery. Taking the total background cross section to be 2 fb from Table [l], we 

show in Table ^ the expected maximum stop discovery mass reach at Tevatron Run II for various amounts 

7 If x? is not a pure bino, its branching ratio to 7G will typically be reduced somewhat when its mass is large (see Fig. 
This will lead to a reduction of the signal cross section at large m^o by typically a few tens of percent. 
8 The diphoton trigger efficiency is close to 100%, so we neglect it here. 

9 Note, however, that the non-negligible mass of the LSP of about 35 GeV assumed in Ref. |3(]] leads to kinematics that 
differ significantly from those in our analysis, in which the gravitino is essentially massless. 
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Figure 3: Cross sections in fb for stop pair production in Run II with t — > cjG, after cuts. No efficiencies 
have yet been applied. The black area is excluded by non-observation of j $t events in Run I. 
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fc 


B 


S for a 5a discovery 


a s x 


Maximum stop mass reach 


2 fb" 1 


4 


14 


7.0 fb 


265 GeV 


4 fb" 1 


8 


18 


4.5 fb 


285 GeV 


15 fb" 1 


30 


31 


2.1 fb 


310 GeV 


30 fb- 1 


60 


42 


1.4 fb 


325 GeV 



Table 2: Number of signal events (S) required for a 5a stop discovery at Tevatron Run II in the j'j'77 $t 
channel and the corresponding signal cross section after cuts and efficiencies and maximum stop mass 
reach. We assume a photon identification efficiency of e 7 = 0.80. The number of background events (B) is 
based on a background cross section of 2 fb from Table ED. 

of integrated luminosity.^ In particular, with 4 fb -1 a stop discovery can be expected in this channel if 
mi < 285 GeV, with S/B of more than 2/1. [^] Including the effects of mixing in the composition of the 
lightest neutralino, a 50% reduction in the Xi ~^ lG branching ratio compared to the pure bino case would 
reduce the stop mass reach by only about 20 GeV. For such a reduction to occur in the relevant neutralino 
mass range of about 200 — 250 GeV, the lightest neutralino would have to be less than half bino. 

5.2 Three-body decay: t -> bW + ^G 

For a large enough splitting between the stop and neutralino masses, the signature of stop pair production 
would beQ jjWW'yy In this analysis, we consider only the dominant hadronic decay mode of both 
W bosons. 

This decay mode of the stop proceeds via three Feynman diagrams, involving an intermediate off-shell 
top quark, chargino or sbottom. We use the full decay matrix elements as given in Ref. |]31|| . Although this 
introduces several additional parameters into the analysis, a few simplifying assumptions may be made 
without becoming too model dependent. For example, assuming that the lightest stop is predominantly 
the superpartner of the right-handed top-quark (£r), eliminates the sbottom exchange diagram altogether. 
Even if the stop contains a mixture of Ir and t^, under our assumption that the lighter stop is the next-to- 
lightest Standard Model superpartner, the sbottom exchange diagram will be suppressed by the necessarily 
10 The cross sections and numbers of events required for discovery are quoted in terms of integrated luminosities at a single 

detector. If data from the CDF and D0 detectors are combined, the integrated luminosity of the machine is effectively doubled. 
11 For comparison, in the case of minimal supergravity a reach of < 180 GeV can be expected in the i — > cxi channel 

with 4 fb at Tevatron Run II; the same reach is obtained in low-energy SUSY breaking scenarios in which the stop is the 

NLSP rather than the neutralino [^l]. In both of these cases, the signal consists of jj $t, with no photons in the final state. 
The backgrounds are again small enough after cuts in this channel that we do not need to tag the b quarks. In the case of 

a discovery, one could imagine tagging the b quarks and reconstructing the W bosons in order to help identify the discovered 

particle. 



15 



larger sbottom mass. As for the chargino exchange, the wino component does not contribute for a tn decay. 
Thus the chargino contribution is dominated by its Higgsino component. Furthermore, if we concentrate 
on scenarios with large values of the supersymmetric mass parameter \i and of the wino mass parameter 
(in which case the charginos are heavy and the neutralino is almost a pure bino), the chargino exchange 
contribution is also suppressed and the dominant decay mode is via the diagram involving an off-shell top 
quark. To simplify our numerical calculations, we have then assumed that only this diagram contributes to 
the stop decay matrix element. We have checked for a few representative points though, that the inclusion 
of the sbottom and chargino diagrams does not significantly change the signal efficiency after cuts as long 
as we require that mg, > m t ~. 

As the W bosons themselves decay, it might be argued that their polarization information needs to 
be retained. However, since we do not consider angular correlations between the decay products, this is 
not a crucial issue; the loss of such information at intermediate steps in the decay does not lead to a 
significant change in the signal efficiency after cuts. This is particularly true for the hadronic decay modes 
of the W, for which the profusion of jets frequently leads to jet overlap, thereby obscuring detailed angular 
correlations. It is thus safe to make the approximation of neglecting the polarization of the W bosons in 
their decay distributions, and we do so in our analysis. 

Before we discuss the signal profile and the backgrounds, let us elaborate on the aforementioned jet 
overlapping. With six quarks in the final state, some of the resultant jets will very often be too close to 
each other to be recognizable as coming from different partons. We simulate this as follows. We count a 
final-state parton (quark or gluon) as a jet only if it has a minimum energy of 5 GeV and lies within the 
pseudorapidity range —3 < rj < 3. We then merge any two jets that fall within a 5R separation of 0.5; the 
momentum of the resultant jet is the sum of the two momenta. We repeat this process iteratively, starting 
with the hardest jet. Our selection cuts are then applied to the (merged) jets that survive this algorithm. 

The signal thus consists of: 

n jets + 27+ $ T (n < 6) 

Hence, all of the SM processes discussed in the previous section yield backgrounds to this signal when up 
to four additional jets are radiated. Now, the radiation of each hard and well separated jet suppresses the 
cross section by a factor of order a s ~ 0.118. Then, since the jj'77 fir backgrounds are already quite 
small after the cuts applied in the previous section (see Table |l]), the backgrounds with additional jets are 
expected to be still smaller. 

There exists a potential exception to the last assertion, namely the background due to ££77 production. 
To get an order of magnitude estimate, the cross section for ti production at Tevatron Run II is 8 pb |32| . 
If both of the photons are required to be energetic and isolated, we would expect a suppression by a factor 
13 For the backgrounds in which one or both of the photons are faked by misidentified jets, we have taken into account the 
larger combinatoric factor that arises when more jets are present to be misidentified. 
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of order a^ m , leading to a cross section of the order of 0.5 fb. This cross section is large enough that we 
need to consider it carefully. 

If both W bosons were to decay hadronically, then a sufficiently large missing transverse energy can 
only come from a mismeasurement of the jet or photon energies. As we have seen in the previous section, 
this missing energy is normally too small to pass our cuts, thereby suppressing the background. If one of 
the W bosons decays leptonically, however, it will yield a sizable amount of ft?. This background can be 
largely eliminated by requiring that no lepton (e or fi) is seen in the detector. This effectively eliminates 
the W decays into e or /i or decays to r followed by leptonic r decays; the remaining background with 
hadronic r decays is naturally quite small without requiring additional cuts. Considerations such as these 
lead us to an appropriate choice of criteria for an event to be selected: 

1. At least four jets, each with a minimum transverse momentum prj > 20 GeV and contained in the 
pseudorapidity interval of —3 < rjj < 3. Any two jets must be separated by SRjj > 0.7. As most of 
the signal events do end up with 4 or more energetic jets (the hardest jets coming typically from the 
W boson decays), this does not cost us in terms of the signal, while reducing the QCD background 
significantly. In addition, the M77 events with both Ws decaying leptonically are reduced to a level 
of order 10 -4 fb by this requirement alone. 

2. Two photons, each with px^ > 20 GeV and pseudorapidity —2.5 < ?? 7 < 2.5. The two photons must 
be separated by at least <5i? 77 > 0.3. 

3. Any photon-jet pair must have a minimum separation of 5Rj~ > 0.5. 

4. A minimum missing transverse energy ftp > 30 GeV. 

5. The event should not contain any isolated lepton with pj- > 10 GeV and lying within the pseudora- 
pidity range — 3 < rj < 3. 

As in the previous section, the cut on ffip serves to eliminate most of the background events with only 
a fake missing transverse momentum (arising out of mismeasurement of jet energies). In association with 
the lepton veto, it also eliminates the bulk of events in which one of the W bosons decays leptonically 
(including the r channel). A perusal of Table ||, which summarizes the major backgrounds after cuts, 
convinces us that the backgrounds to this channel are very small, in fact much smaller than those for the 
previous channel. 

The signal cross section after cuts (but before photon identification efficiencies) is shown in Fig. ||| as a 
function of the stop and Xi masses. We assume that the branching ratio of t — ► bWxi dominates in the 
region of parameter space under consideration. The branching ratio of Xi ~ > lG is again taken from Fig. || 
assuming that X? is a pure bino. The signal efficiency after cuts is about 45%. For small neutralino masses 
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Background 


Cross section after cuts 


after 7 ID 


{jjllZ, Z -> vv) + 2j 


~ 0.003 fb 


~ 0.002 fb 


jjvf/'yy + 2j 


~ 3 x 10" 5 fb 


~ 2 x 10" 5 fb 


(6677, CC77) + 2j 


~ 0.001 fb 


~ 0.0006 fb 


jjn + 2j 


< 0.003 fb 


< 0.002 fb 


firry, WW -> jjjj 


< 10- 4 fb 


< 10~ 4 fb 


M77, VKW — >■ jjiv, t = e, //, or r — ► £ 


~ 0.001 fb 


~ 0.0006 fb 


££77, P^iy — ► 7 jri/, t — > j 


< 0.01 fb 


< 0.006 fb 


Backgrounds with fake photons: 






jj(ee -> 77) + 2j 


~ 7 x 10~ 6 fb 


~ 7 x 10~ 6 fb 


Ji7(j 7) + 2j 


~ 0.02 fb 


~ 0.02 fb 


iiO'j -> 11) + 2j 


~ 0.03 fb 


~ 0.03 fb 


Total 


< 0.07 fb 


~ 0.06 fb 



Table 3: Backgrounds to tt* — > jjWWjj $t with both W bosons decaying hadronically. The photon 
identification efficiency is taken to be e 7 = 0.8 for each real photon. See text for details. 

(m^o <i 50 GeV), though, both the photons and the gravitinos ($r) tend to be soft, leading to a decrease 
in the signal efficiency. For small stop masses (as well as for large stop masses when the stop-neutralino 
mass difference is small), on the other hand, the jets are soft leading to a suppression of the signal cross 
section after cuts. As one would expect, both of these effects are particularly pronounced in the contours 
corresponding to large values of the cross section. The additional distortion of the contours for large 
neutralino masses can, once again, be traced to the suppression of the Xi ~^ lG branching fraction (see 
Fig. §. 

The non-observation of jjWWjj $t events at Run I of the Tevatron already excludes the region of 
parameter space shown in black in Fig. ||. As in the previous section, in this excluded region at least 3 
signal events would have been produced after cuts and detector efficiencies in the 100 pb~ x of Run I data, 
with negligible background. In particular, Run I data excludes stop masses below about 200 GeV, for 
neutralino masses larger than about 50 GeV. 

We show in Table || the expected maximum stop discovery mass reach at Tevatron Run II for various 
amounts of integrated luminosity.^ In particular, with 4 fb -1 , a stop discovery can be expected in this 
channel if m t - < 320 GeVj^] If the lightest neutralino is not a pure bino, the reach at large neutralino 
14 Again, if data from the CDF and D0 detectors are combined, the integrated luminosity of the machine is effectively 
doubled. 

15 For comparison, in the case of minimal supergravity a reach of mj < 190 GeV can be expected in the t — > bWxl channel 
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Figure 4: Cross section in fb for stop pair production with t — ► bW^fG, after cuts. Both W bosons are 
assumed to decay hadronically. The black area is excluded by non-observation of jjWWjj J£t events in 
Run I. 
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fc 


B 


S for a 5(7 discovery 


a s x 4 


Maximum stop mass reach 


2 fb" 1 


0.1 


5 


2.5 fb 


300 GeV 


4 fb" 1 


0.2 


6 


1.5 fb 


320 GeV 


15 fb" 1 


0.9 


8 


0.53 fb 


355 GeV 


30 fb- 1 


1.8 


10 


0.33 fb 


375 GeV 



Table 4: Number of signal events (S) required for a 5a stop discovery at Tevatron Run II in the jjWW'yy / 
Et channel and the corresponding signal cross section after cuts and efficiencies and maximum stop mass 
reach. We take e 7 = 0.80. The number of background events (B) is based on a background cross section 
of 0.06 fb from Table |. 

masses will be reduced. However, the maximum stop mass reach quoted here will not be affected, because 
it occurs for m^o ~ 50 — 100 GeV; in this mass range the lightest neutralino decays virtually 100% of 
the time to due to the kinematic suppression of all other possible decay modes, unless its photino 
component is fine-tuned to be tiny. 

Including a separate analysis of stop production and decay with one or more of the W bosons decaying 
leptonically would yield an increase in the overall signal statistics; however, we do not expect this increase 
to dramatically alter the stop discovery reach. 

5.3 Stop production in top quark decays 

If m ( -+m^o < rrit, then stops can be produced in the decays of top quarks. As we will explain here, most of 
the parameter space in this region is excluded by the non-observation of stop events via direct production 
or in top quark decays at Run I of the Tevatron. However, some interesting parameter space for this decay 
remains allowed after Run I, especially if the lighter stop is predominantly ti,. 

For mi < nib + mw + m x°i so that the stop decays via t — > cxi, the region in which t — > txi is 
possible is almost entirely excluded by the limit on stop pair production at Run I, as shown in Fig. ||. A 
sliver of parameter space in which the stop-neutralino mass splitting is smaller than about 10 GeV remains 
unexcluded. For m 4 - > + mw + m^o, so that the stop decays via i — > bWxi, the signal efficiency in the 
search for direct stop production is degraded for light neutralinos with masses below about 50 GeV and 
for stops lighter than about 150 GeV. This prevents Run I from being sensitive to stop pair production in 
the region of parameter space in which top quark decays to stops are possible with the stops decaying to 
bWxi, as shown in Fig. [|]. In what follows, we focus on this latter region of parameter space. 

As discussed before, if the lightest neutralino is mostly bino, the constraints on its mass are model- 
dependent. The constraints from Tevatron Run I are based on inclusive chargino and neutralino produc- 
with 4 fb" 1 at Tevatron Run II Q. 
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tion ||,|3l| under the assumption of gaugino mass unification; the cross section is dominated by production 
of Xi an d xfx®- If t ne assumption of gaugino mass unification is relaxed, then Run I puts no constraint 
on the mass of Xv At LEP, while the pair production of a pure bino leads to an easily detectable diphoton 
signal, it proceeds only via i-channel selectron exchange. The mass bound on a bino Xi from LEP thus 



depends on the selectron mass [13|. In particular, for selectrons heavier than about 600 GeV, bino masses 
down to 20 GeV are still allowed by the LEP data. If Xi contains a Higgsino admixture, it couples to the 
Z and can be pair-produced at LEP via Z exchange. For an NLSP with mass between 20 and 45 GeV, 
as will be relevant in our top quark decay analysis, the LEP search results limit the Hi component to be 
less than 1%. Such a small H2 admixture has no appreciable effect on the top quark partial width to txi- 
We thus compute the partial width for t — > txi assuming that the neutralino is a pure bino. Taking the 
lighter stop to be t\ = ti cos % + tR sin 9^, we find, 



r(t -» tiB) 



- sin 2 6f H cos 2 6+ 

9 * 36 * 



a Ef 



cos 2 9w m t * B 



where Eg = (m 2 + m 2 ^ — m~)/2mt- The numerical factors in the square brackets in Eq. 5.1 come from 
the hypercharge quantum numbers of the two stop electroweak eigenstates. Clearly, the partial width is 
maximized if the lighter stop is a pure tR state; it drops by a factor of 16 if the lighter stop is a pure ti 
state. In any case, the branching ratio for t — > tB does not exceed 6% for mi > 100 GeV and m^o > 20 
GeVj^ 

The signal from top quark pair production followed by one top quark decaying as in the Standard 
Model and the other decaying to txi, followed by the stop 3-body decay and the neutralino decays to 7G, 
is bbWW'yy $t- This signal is the same (up to kinematics) as that from stop pair production in the 3-body 
decay region. As in the case of stop pair production followed by the 3-body decay, we expect that the 
background to this process can be reduced to a negligible level. Run I can then place 95% confidence level 
exclusion limits on the regions of parameter space in which 3 or more signal events are expected after cuts 



and efficiencies are taken into account. Using the Run I top quark pair production cross section of 6 pb [32| 
and a total luminosity of 100 pb" 1 , we compute the number of signal events as a function of the stop and 
neutralino masses and the stop composition, assuming various values of the signal efficiency after cuts and 
detector efficiencies. If t\ = tR, then Run I excludes most of the parameter space below the kinematic 
limit for this decay even for fairly low signal efficiency ~ 20%, as shown in Fig. ||. If t\ = t^, on the other 
hand, the signal cross section is much smaller and Run I gives no exclusion unless the signal efficiency is 
larger than 75%, which would already be unfeasible including only the identification efficiencies for the two 
photons; even 100% signal efficiency would only yield an exclusion up to ~ 118 GeV. 



At Run II, the top quark pair production cross section is 8 pb [32| and the expected total luminosity 



16 For neutralino masses below mz, as are relevant here, the branching ratio into 7G is virtually 100% (see Fig. g), almost 
independent of the neutralino composition. 
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Stop mass (GeV) 

Figure 5: 95% confidence level exclusion limits for top quark decays to stops from Run I for various signal 
efficiencies, in the case that t\ = tn, which gives the largest event rates. The area to the left of the curves 
is excluded. The case t\ = ti gives no exclusion for the signal efficiencies considered and is not shown here. 
The solid line from upper left to lower right is the kinematic limit for m< = 175 GeV. The solid line from 
upper right to lower left separates the regions in which the 2-body FC decay and the 3-body decay of the 
stop dominate. 
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Figure 6: 5<r discovery contours for top quark decays to stops at Run II with 4 fb _1 for various signal 
efficiencies, in the case that t\ = t^, which gives the smallest event rates. The case t± = tn would be 
discovered virtually up to the kinematic limit even with 10% signal efficiency, and is not shown here. The 
solid lines are as in Fig. ||. 

is considerably higher. This allows top quark decays to txi to be detected for stop and neutralino masses 
above the Run I bound. For t% = £r, top quark decays to stops will be probed virtually up to the kinematic 
limit, even with low signal efficiency ~10% and only 2 fb _1 of integrated luminosity. If i\ = t^, so that 
the signal event rate is minimized, top quark decays to stops would be discovered up to within 10 GeV of 
the kinematic limit for signal efficiencies > 20% and 4 fb _1 of integrated luminosity (see Fig. ^). In this 
region of parameter space, stops would also be discovered in Run II with less than 2 fb _1 via direct stop 
pair production (see Fig. ||). 
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6 Ten degenerate squarks 



Having concentrated until now on the production and decay of the top squark, let us now consider the 
other squarks. In the most general MSSM, the spectrum is of course quite arbitrary. However, low energy 
constraints [34 from flavor changing neutral current processes demand that such squarks be nearly mass 



degenerate, at least those of the same chirality. Interestingly, in many theoretical scenarios, such as the 
minimal gauge mediated models, this mass degeneracy between squarks of the same chirality happens 
naturally; in addition the mass splitting between the left-handed and right-handed squarks associated with 
the five light quarks turns out to be small. For simplicity, then, we will work under the approximation 
that all of these 10 squarks (namely, ul,r, d,L t n, &l,r, sl,r and 1)l,r) are exactly degenerate. 

6.1 Production at Tevatron Run II 

While the cross sections for the individual pair-production of the cl,r, sl,r and ol,r are essentially the 
same as that for a top-squark of the same mass, the situation is more complicated for squarks of the 
first generation. The latter depend sensitively on the gluino mass because of the presence of i-channel 
diagrams. Moreover, processes such as uu — > ulu r or dd — > di t ndL,R become possible and relevant. Of 
course, in the limit of very large gluino mass, the squark production processes are driven essentially by QCD 
and dominated by the production of pairs of mass eigenstates, analogous to the top squark production 
considered already. In particular, at the leading order, the total production cross section for the ten 
degenerate squarks of a given mass is simply ten times the corresponding top squark production cross 
section. 

Since a relatively light gluino only serves to increase the total cross section (see Fig. |7|) , it can be argued 
that the heavy gluino limit is a conservative one. To avoid considering an additional free parameter, we 
shall perform our analysis in this limit. To a first approximation, the signal cross sections presented below 
will scalej^] with the gluino mass approximately as shown in Fig. |7[ 

Like top squarks, the 10 degenerate squarks can also be produced via cascade decays of heavier super- 
symmetric particles. To be conservative, we again neglect this source of squark production by assuming 
that the masses of the heavier supersymmetric particles are large enough that their production rate at 
Tevatron energies can be neglected. 

The NLO cross sections for production of ten degenerate squarks including QCD and SUSY-QCD 
corrections have been implemented numerically in PROSPINO |18[| . We generate squark production events 
using the LO cross section evaluated at the scale \x = nig, improved by the NLO K-factor obtained from 
17 That the gluino exchange diagram has a different topology as compared to the (dominant) quark-initiated QCD diagram 
indicates that corresponding angular distributions would be somewhat different. Thus, the efficiency after cuts is not expected 
to be strictly independent of the gluino mass. For the most part, though, this is only a subleading effect. 
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1000 10000 
Gluino mass (GeV) 

Figure 7: Gluino mass dependence of the NLO cross section for production of ten degenerate squarks in 
2 TeV pp collisions, from PROSPINO pl| . Shown are the cross sections for qq* production normalized to 
the value at large m g , for common squark masses of 200, 300, 400 and 500 GeV. Production of qq is small 
at a pp collider and is neglected here; it yields an additional 3-15% increase in the total cross section at 
low m g for this range of squark masses. Cross sections are evaluated at the scale [i = rriq. 
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Figure 8: LO and NLO cross section for the production of ten degenerate squarks in the heavy gluino limit 



in pp collisions at Tevatron Run I (1.8 TeV) and Run II (2.0 TeV), from PROSPINO [18|. Cross sections 



are evaluated at the scale \x 



PROSPINO |18| (see Fig. ||), in the limit that the gluino is very heavy. The K-factor varies between 1 and 
1.25 for niq decreasing from 550 to 200 GeV. As in the case of the top squark analysis, we use the CTEQ5 
parton distribution functions [19] and neglect the shift in the pt distribution of the squarks due to gluon 
radiation at NLO. 



6.2 Signals in low-energy SUSY breaking 

The decays of the ten degenerate squarks are very simple. As long as they are heavier than the lightest 
neutralinoj^ they decay via q — > qx± — ► <TfG. The signal and backgrounds are then identical to those of 



the two-body FC stop decay discussed in Sec. [5JJ , and consequently we use the same selection cuts. In 
fact, in view of the tenfold increase in the signal strength, we could afford more stringent cuts so as to 



18 As before, we do not allow the possibility of cascade decays through other neutralinos/charginos. Were we to allow these, 
the channel we are considering would be somewhat suppressed, but additional, more spectacular, channels would open up. 
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eliminate virtually all backgrounds, but this is not quite necessary. 

The signal cross section after cuts (but before efficiencies) for production of ten degenerate squarks in 
the heavy gluino limit is shown in Fig. ^ as contours in the rriq-m^o plane. We assume that the squark q 
decays predominantly into qxi- The branching ratio of Xi ~^ lG ls taken from Fig. ^ assuming that Xi 
is a pure bino. Clearly, the effect of the kinematic cuts on the signal is very similar to that in the case of 
the 2-body FC decay of the stop. The mass reach, of course, is much larger due to the tenfold increase in 
the total cross section; also, unlike in the case of the stop, the 2-body decay is dominant throughout the 
entire parameter space. 

The non-observation of jj'77 J£t events at Run I of the Tevatron excludes the region of parameter space 
shown in black in Fig. ^. As in our stop analysis, in this excluded region, at least 3 signal events would 
have been detected in the 100 pb _1 of Run I data, with negligible background. In particular, we estimate 
that Run I data excludes the ten degenerate squarks up to a common mass of about 280 GeV. As in the 
case of 2-body FC stop decays, when the mass splitting between the squarks and the neutralino is too small 
(i.e., less than about 40 GeV), the jets become soft and the signal efficiency decreases dramatically, leaving 
an unexcluded region of parameter space. The D0 search for inclusive pp — » x% + X with X2 ~^ 7Xi HI 
yields a limit on the production cross section of about 0.5 pb for parent squark masses above about 250 
GeV. Interpreted in terms of pair production of 10 degenerate squarks in the large m g limit with X2 ~^ 1X\ 
reidentified as Xi ~^ 7^ again increases the signal efficiency by a factor of ~ 2.7 because every event 
contains two photons |3(|; the D0 analysis then yields a bound on the common squark mass of about 275 
GeV, again in rough agreement with our result. PI 

Taking the total background cross section to be 2 fb (see Table [l]), we show in Table [5] the expected 
maximum discovery mass reach for ten degenerate squarks at Tevatron Run II for various amounts of 
integrated luminosity.0 In particular, with 4 fb -1 a squark discovery can be expected in this channel if 
mg < 360 GeV. Again, if the lightest neutralino is not a pure bino, the reach at large neutralino masses 
will be reduced. However, this will have very little effect on the maximum squark mass reach quoted here, 
because this maximum reach occurs for m^o ~ 100 — 150 GeV, where all neutralino decay modes other 
than 7G suffer a large kinematic suppression. From Fig. B it is evident that the neutralino branching 
fraction into 7G will be reduced by no more than 10% in this mass range as long as the neutralino is at 
least 50% bino; such a reduction in the neutralino branching fraction will lead to a reduction of only a few 
GeV in the maximum squark mass reach. 



19 Ref. fx) quotes a squark mass bound of 320 GeV for the case rriq <C m g in the context of low-energy SUSY breaking; 
we expect that this is due to the contribution of chargino and neutralino production to the total SUSY cross section in their 
analysis. 

20 Again, if data from the CDF and D0 detectors are combined, the integrated luminosity of the machine is effectively 
doubled. 
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Figure 9: Cross section in fb for production of 10 degenerate squarks in Run II in the heavy gluino limit 
with q — ► q'jG, after cuts. The black area is excluded by non-observation of jj'77 $t events in Run I. 
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fc 


B 


S for a 5a discovery 


a s x 


Maximum squark mass reach 


2 fb^ 1 


4 


14 


7.0 fb 


345 GeV 


4 fb" 1 


8 


18 


4.5 fb 


360 GeV 


15 ft" 1 


30 


31 


2.1 fb 


390 GeV 


30 fb- 1 


60 


42 


1.4 fb 


405 GeV 



Table 5: Number of signal events (S) required for a 5a squark discovery at Tevatron Run II, assuming 
production of 10 degenerate squarks in the limit that the gluino is very heavy, and the corresponding signal 
cross section after cuts and efficiencies and maximum squark mass reach. We take e 7 = 0.80. The number 
of background events (B) is based on a background cross section of 2 fb from Table 0. 

7 Conclusions 

In models of low-energy SUSY breaking, signatures of SUSY particle production generically contain two 
hard photons plus missing energy due to the decays of the two neutralino NLSPs produced in the decay 
chains. Standard Model backgrounds to such signals are naturally small at Run II of the Tevatron. We 
studied the production and decay of top squarks at the Tevatron in such models in the case where the 
lightest Standard Model superpartner is a light neutralino that predominantly decays into a photon and a 
light gravitino. We considered 2-body flavor-changing and 3-body decays of the top squarks. The reach of 
the Tevatron in such models is larger than in the standard supergravity models and than in models with 
low-energy SUSY breaking in which the stop is the NLSP, rather than the neutralino. We estimate that 
top squarks with masses below about 200 GeV can be excluded based on Run I data, assuming that 50 
GeV < m^o < mi — 10 GeV. For a modest final Run II luminosity of 4 fb -1 , stop masses up to 285 GeV 
are accessible in the 2-body decay mode, and up to 320 GeV in the 3-body decay mode. 

Top squarks can also be produced in top quark decays. We found that, within the context of low-energy 
SUSY breaking with the stop as the next-to-next-to-lightest SUSY particle, the region of parameter space 
in which stop production in top quark decays is possible is almost entirely excluded by Run I data if the 
lighter stop is predominantly right handed; however, an interesting region is still allowed if the lighter stop 
is predominantly left handed, due to the smaller branching ratio of t — ► t^Xi- Run II will cover the entire 
parameter space in which top decays to stop are possible. 

We also studied the production and decay of the ten squarks associated with the five light quarks, 
assumed to be degenerate. In models of low-energy SUSY breaking, the decays of the ten degenerate 
squarks lead to signals identical to those for the 2-body flavor-changing stop decays. The cross section for 
production of ten degenerate squarks at the Tevatron is significantly larger than that of the top squark. 
We estimate that the 10 degenerate squarks with masses below about 280 GeV can be excluded based on 
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Run I data, assuming that m^o < rriq — 40 GeV. For a final Run II luminosity of 4 fb , squark masses 
as large as 360 GeV are easily accessible in the limit that the gluino is very heavy. The production cross 
section, and hence the discovery reach, increases further with decreasing gluino mass. 
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